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Abstract 


An approach for an optimization-based integrated controls- structures 
design is presented for a class of flexible spacecraft that require fine 
attitude pointing and vibration suppression . The integrated design problem 
is posed in the form of simultaneous optimization of both structural and 
control design variables . The approach is demonstrated by application to 
the integrated design of a generic space platform and to a model of a 
ground-based flexible structure. The numerical results obtained indicate 
that the integrated design approach can yield spacecraft designs that have 
substantially superior performance over a conventional design wherein the 
structural and control designs are performed sequentially. For example , 
a 40 -percent reduction in the pointing error is observed along with a 
slight reduction in mass, or an almost twofold increase in the controlled 
performance is indicated with more than a 5-percent reduction in the 
overall mass of the spacecraft (a reduction of hundreds of kilograms). 


Introduction 

Current spacecraft design customarily involves an 
iterative sequence of designs performed separately 
within the structural and control disciplines. Struc- 
tural design is performed first and is based on loading 
considerations such as would occur during launch, re- 
boost, or component operational maneuvers. Sizes 
and masses of mission-related components are esti- 
mated and a structure that maintains the desired 
component relationships during operations is de- 
signed. A controller is next designed for the fixed 
structure to orient, guide, and/or move the space- 
craft to obtain the required performance. The con- 
trol design must also provide satisfactory closed-loop 
stability and robustness properties. If the nominal 
structure does not admit an adequate control de- 
sign, it is returned to the structural discipline for 
modification. After modification, the structure is re- 
turned to the control discipline for redesign. This 
iterative process continues until a satisfactory com- 
promise is found between the mission and control re- 
quirements. This separate-discipline approach has 
been successfully used in the past and works well 
in cases where a relatively high-stiffness structure is 
attainable and where nonstructural components are 
concentrated masses and inertias, or where perfor- 
mance requirements are not stringent. However, fu- 
ture space structures and performance requirements 
do not fit this category. 

A number of future space missions will use large 
flexible structures in low-Earth and geostationary or- 
bits. Possible structures include space science plat- 
forms, space processing facilities, and Earth observa- 
tion systems. Such structures typically require large 
distributed-mass components such as booms, solar 


arrays, and antennas whose dimensions range from a 
few meters to possibly hundreds of meters. To mini- 
mize the costs of construction, launching, and opera- 
tions, it is necessary to make the structure as light as 
possible. However, the combination of large size and 
low structural mass leads to increased flexibility and 
makes more difficult the control of the structure and 
its components to a specified precision in attitude 
and shape. 

The combination of flexibility and low mass cre- 
ates a special dynamical interaction problem that is 
characteristic of large flexible space structures. In 
the design of controllers for physical systems, some 
trade-off almost always has to be performed between 
design-model accuracy and mathematical complex- 
ity. The more accurate analysis models often re- 
quire computational time that is too excessive to 
qualify them as design models for control purposes. 
Also, many of the most widely-used multivariable 
controller design techniques work best for moderate- 
order, linear, time- invariant design models. Typi- 
cally, in practice, high-order nonlinear models are 
linearized about some operating condition and have 
their model order reduced to produce design models 
that conform to computational limitations or com- 
pensator implementation constraints. These prac- 
ticalities introduce modeling errors in the form of 
unmodeled dynamics that must be accounted for in 
the controller design process. Space structure con- 
troller design models are generally found through 
some order-reduction procedure applied to a high- 
order analysis model obtained from finite element 
techniques. The order-reduction process essentially 
deletes a portion of the finite element model to pro- 
duce a lower-order controller design model. Although 


no longer contained in the design model, the un- 
modeled dynamics, represented by the deleted por- 
tion, still appear in the measurement vector and 
therefore can be affected by model-based control in- 
puts. Care must be taken in the design of controllers 
for lightly damped systems that the control inputs 
do not spillover into the unmodeled dynamics with 
disastrous destabilizing effects (ref. 1). 

Controls-structures interaction (CSI) in the form 
of destabilizing spillover has been verified in sim- 
ple Earth-based laboratory experiments as well as in 
the design, analysis, ground development and test, 
and/or flight operation of space systems in industry 
(ref. 2). The current approach to solving CSI prob- 
lems is to design the spacecraft to avoid undesired 
dynamical interaction. This approach generally re- 
quires either stiffening the structure or slowing down 
the control system response. Stiffening the structure 
simplifies the control design problem in that the pre- 
dominant dynamics tend toward rigid body but is 
costly in terms of mass as well as launch packaging 
and weight, leading to increased fuel consumption. 
Slowing down the control response produces control 
inputs that have less of a chance of producing de- 
stabilizing effects but is costly in terms of reduced 
performance capability. Neither approach is com- 
pletely satisfactory. What is needed is a new T design 
approach that is capable of avoiding the damaging as- 
pects of CSI while identifying and taking advantage 
of any beneficial aspects. 

The development and experimental validation of 
such an approach is a primary purpose of NASA’s 
CSI program (ref. 3). Fundamentally, the CSI pro- 
gram recognizes the high degree of coupling that ex- 
ists between the control and structural disciplines 
when dealing with flexible structures. For example, 
controllers designed to be robust to unmodeled dy- 
namics may need very low T gain (and therefore, pro- 
duce low performance) (ref. 4). In order to main- 
tain high performance, it would then be necessary 
to redesign the structure to increase the frequen- 
cies of the higher modes likely to be affected by 
spillover. If the thickness of a structural member 
is changed, the dynamics will change, which would 
then change the control law and the required ac- 
tuator sizes (and masses). These changes would, 
in turn, again change the structural model. One 
of the goals of the CSI program is to integrate the 
structures and controls design disciplines. Rather 
than design control and structural aspects of space 
structures separately in the sequential manner de- 
scribed earlier, approaches are being developed for 
a unified controls-structures modeling, analysis-and- 
design method that would allow a complete iteration 


on all critical (control and structure) design variables 
in a single integrated computational framework. 

This paper considers missions that include flex- 
ible spacecraft with no articulated appendages and 
require fine attitude pointing and vibration suppres- 
sion (e.g., large space antennas). Missions that in- 
volve flexible spacecraft with articulated multiple 
payloads, missions requiring large- angle slewing of 
spacecraft without articulated appendages, and mis- 
sions that include general nonlinear motion of a flexi- 
ble spacecraft with articulated appendages and robot 
arms are not considered. Two controller strategies 
are considered, namely, the static and dynamic dis- 
sipative controllers, which are well known for their 
stability robustness in the presence of unmodeled dy- 
namics and parametric uncertainties (refs. 4 and 5). 
The approach selected for the development of inte- 
grated controls-structures design methodology is an 
optimization-based procedure employing mathemat- 
ical programming techniques. The optimization ap- 
proach allows freedom and variety in selecting the 
potentially large number of design variables. The op- 
timization approach is also the one commonly used 
in the field of structural design. In fact, many papers 
describing optimization-based integrated design have 
appeared in the literature during the last 12 years. 
An excellent literature survey can be found in ref- 
erence 6. However, only recently have applications 
gone beyond analytical studies with fictitious models, 
relatively simple laboratory apparatus, or unrealistic 
control designs. Some of the recent efforts include 
the works by Balakrishnan using infinite-dimensional 
models (ref. 7); by Bossi, Hunziker, and Kraft us- 
ing Q-parameterization theory along with shape op- 
timization (ref. 8); the studies by Maghami et al. 
on the integrated design of large structures (with 
thousands of degrees of freedom) using dissipative 
compensators (refs. 9-11); and a paper by Woodard 
et al. on integrated optimization with varying actu- 
ator mass (ref. 12). One of the unique features of 
the CSI program is that integrated designs will be 
fabricated and the methodology validated in the lab- 
oratory. This feature places the constraint of realism 
on the design process, which is not normally found 
elsewhere. 

The purpose of this paper is to describe a por- 
tion of the integrated design methodology and results 
from the CSI program. This work differs from past 
studies in several ways. The design treats large re- 
alistic applications* a derivative of a geostationary 
platform designed in support of Earth observation 
sciences (ref. 13) and the initial phase of an evolution- 
ary laboratory structure in the CSI program (ref. 14). 

A design similar to the integrated design reported 
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herein and that uses the CSI model is actually being 
constructed at the Langley Research Center (LaRC) 
for the purpose of validating the integrated design 
methodology. Also, we recognize that the effect of 
the control system on the structural mass and con- 
straints can be significant when the magnitude of the 
control effort determines the mass of the control sys- 
tem and its power supply (ref. 15). This realistic ef- 
fect is accounted for by the inclusion of a mathemat- 
ical expression that couples the actuator masses and 
the controller gains. Finally, the controller formu- 
lation uses recently developed theory for dissipative 
controllers not previously employed in applications 
of this size. 

We begin with a description of the mathematical 
model of the structure and controller design methods. 
Next, we offer an overview of the optimization-based 
integrated design procedure employed in this paper. 
Afterwards, descriptions of the structure, design vari- 
ables and constraints, and integrated design problem 
statements are given for each of the applications. Fi- 
nally, results from the integrated design process are 
discussed. An appendix giving the fundamentals of 
the dissipative controller theory is provided. 

Mathematical Model of the Plant 

The linear, time-invariant, mathematical model 
of a flexible space structure is given by 

Jx + Dx + Kx = Bu (1) 


where q r is an r x 1 vector of modal amplitudes; 
Jr, Dr, K r are, respectively, the generalized iner- 
tia, damping, and stiffness matrices; and <F r is an 
r x n matrix whose columns are the r open-loop 
eigenvectors associated with the included modes. If 
the mode shapes are normalized with respect to 
the inertia matrix, and modal damping is assumed, 
then Jr = Irxr , D r = Diag (2(qltq , 2(^2^2, • ■ ’ ? 2£r^r ) , 
and K r = Diag (u > \ , , . . . , } , where uji and 

Q (i = 1, 2, . . . , r) are the open-loop frequencies and 
damping ratios, respectively. Here, the sensors 
are collocated and compatible with the actuators, 
whereby 


y p = r q r 1 
y r = Tq r j 


(4) 


Note that the collocation of the sensors and actuators 
is necessary for the implementation of the dissipative 
controllers. (See the appendix.) This collocation 
guarantees the system to be minimum phase (i.e., 
all transmission zeros are in the left-half plane) and 
will enhance the stability robustness of the overall 
system. 

In defining the state vector z, 


z=(q r q r ) T (5) 


the dynamics of the system, assuming no appreciable 
sensor and/or actuator dynamics, can be written in 
a first-order form 


yp — Cpx 

yr = C r x 


( 2 ) 


where x is an n x 1 displacement vector, J is the 
positive-definite inertia matrix, D is the open-loop 
damping matrix, K is the nonnegative-definite stiff- 
ness matrix, B is an n x m control influence matrix, 
u is an m x 1 control input vector, y v and y r are, 
respectively, ! x 1 position and rate measurements, 
and C p and C r are the corresponding output in- 
fluence matrices. The second-order representation 
of the structure given in equation (1) is obtained 
through finite element modeling. The order of a 
large flexible space structure (LFSS) can be quite 
large. For design and analysis purposes, the order 
of the system is reduced to a design size. This re- 
duction is accomplished by using a modal truncation 
approach wherein every mode in the controller band- 
width is retained and the remaining modes are trun- 
cated. The system equations in modal coordinates 
for the retained modes are written as 


J r q r + D r q r + K r q r = $[Bue T T u (3) 


z = Az + Bu 


(6) 


where 


A = 

Orxr 

-K r 

B = 

V 
> 


and 

y P = [ T 0 ] z = C p z 
y r — [ 0 r] z = C r z 

The general configuration of the controlled structure 
is illustrated in the block diagram shown in figure 1. 
In this figure, w, v p , and v r represent disturbances 
at the input, position output, and the rate output, 
respectively; r denotes the reference signal vector; 
and y per refers to those outputs that are used for 
performance evaluations and not in the feedback 
loop. In this configuration, the system equations are 
written as 

z — Az + Bu -|- Hw (7) 
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y p — C p z + SpVp 

y r = C r z + S r v r (8) 

yper = Cper Z , 

where H, Sp, and S r are influence matrices for the 
disturbances w, v p , and v r , respectively, and C per 
is the output influence matrix associated with y per . 
It is assumed that w, v p , and v r are uncorrelated, 
zero-mean stationary white-noise processes with unit 
intensities. 

Controller Design Methods 

Control system design for LFSS’s is a challenging 
problem because of their special dynamic character- 
istics, which include the following — a large number 
of structural modes wflthin the controller bandwidth; 
low, closely spaced structural mode frequencies; very 
small inherent damping; and lack of accurate knowl- 
edge of the parameters. To implement the con- 
troller, it must be of a reasonably low^ order and must 
also satisfy the performance specifications (i.e., root 
mean square (rms) pointing error and closed-loop 
bandwidth). It must also have robustness to non- 
parametric uncertainties (i.e., unmodeled structural 
modes), and to parametric uncertainties (i.e., errors 
in the knowledge of the design model). Two ma- 
jor categories of controller design methods for LFSS 
are model-based controllers (MBC’s) and dissipative 
controllers. An MBC generally consists of a state 
estimator (a Kalman-Bucy filter or an observer) fol- 
lowed by a linear-quadratic regulator. The state esti- 
mator uses the knowledge of the design model (con- 
sisting of the rotational rigid-body modes and a few 
elastic modes) in its prediction part. Using multi- 
variable frequency-domain design methods, such con- 
trollers can be made robust to unmodeled structural 
dynamics; that is, the spillover effect can be over- 
come (ref. 4). However, such controllers generally 
tend to be very sensitive to uncertainties in the design 
model — in particular, to uncertainty in the structural 
mode frequencies (refs. 4 and 5). An analytical ex- 
planation of this instability mechanism may be found 
in reference 5. Achieving robustness to real paramet- 
ric uncertainties is, as yet, an unsolved problem, al- 
though considerable research activity is in progress 
in that area. 

Static Dissipative Controller 

In view of the sensitivity problem of MBC’s, dissi- 
pative controllers, which use collocated and compat- 
ible actuators and sensors, offer an attractive alter- 
native. Details of dissipative controller theory may 


be found in the appendix. Dissipative controllers 
use special passivity-type input /output properties of 
the plant and offer robust stability in the presence 
of both nonparametric and parametric uncertainties 
(refs. 4 and 5). The simplest controller of this type 
is the constant-gain dissipative controller. Using col- 
located torque actuators and attitude and rate sen- 
sors, the constant-gain dissipative control law is given 

by: j 

u = - G p y p - G r y r (9) ! 

where y p and y r are the measured (3m x 1, where m J 

is the number of three- axis sensors) attitude and rate; } 

Gp and G r are 3m x 3m symmetric, positive-definite 5 

gain matrices. This control law has been proven to i 

give guaranteed closed-loop stability despite unmod- 
eled elastic modes, parameter errors, certain types 
of actuator and sensor nonlinearities (such as satura- 
tion and dead zone; see the appendix for details), and 
first-order actuator dynamics (ref. 4). Robust stabil- 
ity is generally not guaranteed if the actuator dy- 
namics are second or higher order, or in the presence 1 

of digitization, time delays, and control hysteresis. I 

However, if acceleration feedback is permissible, it \ 

is possible to obtain guaranteed stability for second- = 

order actuator dynamics (refs. 5 and 16). For practi- * 

cal space implementation, rate measurements may be I 

obtained from rate gyros and attitude measurements 1 

may be provided by star sensors. However, if linear 
velocity is required (e.g., in the case of ground-based | 

experiments) velocity signals can be obtained by inte- | 

grating accelerometer signals with the aid of washout 1 

filters to asymptotically remove the accelerometer | 

bias. The inclusion of washout filters would, how- 
ever, hinder the dissipativity of this control design. 

Finally, the drawback of this controller is that the 
achievable performance is inherently limited because 
of its simple mathematical structure. 

Dynamic Dissipative Controller 

To obtain higher performance and still retain the 
highly desirable robust stability, dynamic dissipative 
compensators can be used. The main characteris- 
tic of all dissipative controllers is that they do not 
rely on the knowledge of the design model to ensure 
stability, although they use it to obtain the best pos- 
sible performance. An n c -order (two-level) dynamic 
dissipative controller is given by: 

x c = A c x c + B f y r (10) 

u c = — Gx c — G r y r — G p y p (1 1) 

where A C ,B C , and G are the compensator system, 
input, and output matrices, respectively, and G r and 
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Gp are symmetric, positive-definite rate gain and 
position gain matrices. Furthermore, A c is strictly 
Hurwitz, and the Kalman-Yacubovich relations hold 
(ref. 17): 

A<fP + PA C = — Q (12) 

G = B^P (13) 

where P = P T > 0 and Q = Q r > 0. 

Equations (10)— (13) represent a two-level con- 
troller wherein the inner loop consists of static 
position-plus-rate feedback and the outer loop con- 
sists of a dynamic compensator. This controller 
assures robust asymptotic stability regardless of 
unmodeled structural dynamics or parametric un- 
certainties (refs. 5 and 16). In the absence of zero- 
frequency rigid-body modes (e.g., for a ground-based 
experiment), and G r can be null matrices with- 
out destroying the robust asymptotic stability; that 
is, the inner loop is not required. Further results 
for systems with zero- frequency modes are given in 
the appendix. As with the static dissipative con- 
troller, this control law has been proven to give guar- 
anteed closed-loop stability despite unmodeled elas- 
tic modes and parameter errors. However, robust 
stability is generally not guaranteed if the actuator 
dynamics, nonlinearities, and time delays are present 
or if digital implementation is used. Also as with the 
static dissipative controller — for practical space im- 
plementation, rate measurements may be obtained 
from rate gyros and attitude measurements may be 
provided by star sensors. However, if linear veloc- 
ity is required (e.g., in the case of ground-based ex- 
periments), velocity signals can be obtained by inte- 
grating accelerometer signals with the aid of washout 
filters to asymptotically remove the accelerometer 
bias. The inclusion of washout filters would, however, 
hinder the dissipativity of this control design. 

Optimization-Based CSI Design 

The methodology employed in this study uses op- 
timization techniques to integrate the separate disci- 
plines of control and structural design. Fundamen- 
tally, the integrated design goals and requirements 
are formulated as a nonlinear programming prob- 
lem wherein selected control-structure design vari- 
ables are chosen, subject to structural and perfor- 
mance constraints, to optimize an objective function 
whose construction captures the essence of the inte- 
grated design problem. It is felt that this approach 
takes advantage of the highly-developed optimization 
technology in both the controls and structures fields 
and provides the designer an opportunity to obtain 
nonintuitive results. 


Nominally, an initial structure-control configura- 
tion defines a set of design variables that are avail- 
able for modification. Typical controller design vari- 
ables are size and location of sensors and actuators, 
and elements of dynamic compensator system ma- 
trices. Structural design parameters could be the di- 
mensions of structural members along with mass and 
densities of materials. An example of a realistic inte- 
grated design problem would be, for some given input 
or maneuver, to design a spacecraft that provides the 
best possible reductions in dynamic response ampli- 
tude (below the amplitudes of a nominal spacecraft) 
with no increase in structural mass or control system 
energy requirements. As indicated by the preced- 
ing examples, most integrated design objectives and 
constraints can be expected to be dependent func- 
tions of the design variables that are obtained as the 
outputs of some separately developed computer con- 
trol synthesis (ref. 18) or structural analysis pack- 
ages (ref. 19). Integrated design is then a compu- 
tationally intensive process requiring supercomputer 
technology. 

A flowchart summarizing the process for optimi- 
zation-based integrated design is given in figure 2. 
A user- provided executive program drives an opti- 
mization module that, at each stage, (1) uses the 
current values of the design variables to compute the 
objective function and constraints, (2) evaluates the 
objective function and constraints for optimality and 
feasibility, and (3) if necessary, generates a new set 
of design variables defining a structure-control con- 
figuration with better integrated-design characteris- 
tics. Aspects of the computer implementation of this 
process are discussed in reference 20. 

A version of the integrated design software tool 
CSI-DESIGN, which is under development at the 
LaRC, was used to perform the numerical studies re- 
ported herein. CSI-DESIGN is intended for research 
purposes and is composed of public domain software. 
It uses in-core data transfer and is ultimately limited 
in the size of problems it can successfully treat. The 
package has linked control, structural, and optimiza- 
tion modules as shown in figure 2. The control com- 
pensator models used in this study were static and 
dynamic dissipative. The benefits and motivation for 
using this control strategy for the control of flexible 
structures are discussed in the appendix. A descrip- 
tion of the contents of the CSI-DESIGN structural 
module may be found in reference 20. Integrated 
optimization was performed with a four-processor 
Alliant FX/80 digital computer that used the Au- 
tomated Design Synthesis (ADS) software (ref. 21). 
The interior penalty function method of ADS was 
used to solve the nonlinear programming problems. 
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In this method, the constrained optimization prob- 
lem is transformed into an unconstrained problem 
through creation of a pseudo-objective function that 
is the sum of the original objective function and an 
imposed penalty function, which is a function of the 
constraints (ref. 22). The Reverse-Cuthill-McKee al- 
gorithm for minimizing the bandwidth of the banded 
stiffness and mass matrices was used to reduce com- 
putational requirements (ref. 23). Additionally, ana- 
lytical expressions for eigenvalue and eigenvector sen- 
sitivity (with respect to the structural design vari- 
ables) (ref. 24) were used in the integrated design 
process to approximate the eigenvalues and eigenvec- 
tors at design points that are in the neighborhood of 
the nominal design point. This approximation was in 
the form of a first-order Taylor scries approximation 
and resulted in substantial computational savings be- 
cause it removes the need for costly computation of 
structural eigenvalues and eigenvectors at many of 
the optimization moves. 

Integrated Design of the EPS Model 

The Earth pointing system (EPS), a generic 
model of the geostationary platform is a multiuser 
concept and consists of a 10-bay, 30-tn-long truss 
structure with two radial rib antennas (7.5 m and 
15 m diameter) at the ends (fig. 3). All the members 
(i.c., constituting the truss, the antennas, and the an- 
tenna supports) are assumed to be hollow tubes with 
circular cross sections and thicknesses oEI.59 mm. 
The antennas are fixed with respect to the truss, so 
the problem is that of controlling the pointing and 
vibration of the entire structure. It is assumed that 
a three-axis control moment gyro (CMG) and collo- 
cated attitude and rate sensors, located at the center- 
of-mass of the structure, are used for accomplish- 
ing the control. No sensor and actuator dynamics 
are considered. The modal frequencies for the first 
10 flexible modes are given in table 1. The first flex- 
ible mode had a frequency of 0.58 Hz, corresponding 
to a mode of the large antenna support structure, and 
the first flexible mode of the truss structure (mode 7) 
had a frequency of about 6.6 Hz. A design model 
consisting of three rigid-body modes and the first 
10 flexible modes of the structure was used in the 
design process. An open-loop modal damping ratio 
of 0.5 percent was assumed for the flexible modes. 

The integrated design problem w r as formulated as 
a single-objective optimization problem. The struc- 
tural design variables used were outer diameters of 
the truss and antenna support members with the 
thicknesses fixed. In particular, the truss was di- 
vided into three sections shown in figure 3, and the 


outer diameters of the longerons, battens, and diag- 
onals within each section constituted nine structural 
design variables. Two additional structural design 
variables were the outer diameters of the support 
members for the two antennas, thus constituting a 
total of 11 structural design variables. The num- 
ber of control design variables varied and depended 
on the type of controller used- For the static dissi- 
pative design, the elements of the Cholesky factor- 
ization matrices of the positive-definite position and 
rate gain matrices G p and G r are chosen as design 
variables: 

G p = LpLp 1 

T (14) 

G r - L r Lj J 

In dynamic dissipative designs, the elements of the 
compensator matrices A c and B c , the elements 
of weighting matrix Q, and the elements of the 
Cholesky factorization matrices Lp and L r are used 
as design variables. Within this configuration, the 
following two design problems are considered. 

Design Problem I 

The system configuration for the first design prob- 
lem is illustrated in the block diagram showm in fig- 
ure 4. Here, white-noise disturbances of unit inten- 
sity are applied to the structure at the same locations 
as the control inputs, and no disturbances at the out- 
puts are considered. The system equations are given 
as 


z — Az -j- Bu T Bw 

(15) 

y p = C p z \ 


y r = C r z > 

(16) 

Yper — e ~ Gp er z J 



In this design problem, the closed-loop performance 
measure is the steady-state rms pointing error vec- 
tor e at the large antenna due to white-noise distur- 
bance of unit intensity at the inputs. To achieve a 
realistic design, constraints are placed on the steady- 
state average control power and the total structural 
mass. Additional side constraints are placed on the 
structural design variables for safety and practicality 
reasons. Lower bounds are placed on these variables 
to satisfy structural integrity requirements against 
buckling and stress failures. On the other hand, up- 
per bounds are placed on these variables to accom- 
modate manufacturing limitations. The first design 
problem is described as follows. 

Minimize the steady-state rms pointing error at 
the large antenna 

J^lim (TV {£[e(i)e r (*)]}) 1/2 (17) 
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with respect to the tube outer diameters d,{i = 
1, 2, . . . , 11), and the elements of L p and L r (where 
L p and L r are 3x3 lower- triangular Cholesky factors 
of G p and G r given in eq. (14)), subject to the 
constraints: 


lim (Tr {# [u (<) u T (<)] }) < Tmax (18) 

A/ S tr < A/max 


d m in — di < d m ax (^0) 

where e is the 3x1 attitude vector at the large 
antenna, Tr( ) denotes trace of (),£() represents 
expected value of ( ), and Af s tr denotes the structural 
mass. The subscripts min and max define the range 
of allowable values. 

The steady-state covariance of the pointing er- 
ror vector P e is computed from the steady-state 
covariance of the state vector; that is, 

P e ee hm Tr [e (ee T )] = Tr (c per PzCj er ) (21) 
where P z = lim E( zz r ) denotes the steady-state 

t—+o O > 

covariance of the state that is determined from the 
solution of the following Liapunov equation (ref. 25): 

A C ]P Z + P Z A^ = -BB r (22) 

in which A c i is the closed-loop state matrix and B is 
the control (input) influence matrix. 

The results for this design problem are summa- 
rized in table 2. An initial control design was first ob- 
tained for the nominal structure. Using a pole place- 
ment technique, a decoupled controller was designed 
to provide adequate performance for each of the three 
rigid-body modes, with no consideration given to the 
flexible modes. That is, the static-dissipative control 
gain matrices were diagonal, with elements chosen 
to give satisfactory closed-loop frequency and damp- 
ing for the rigid-body dynamics and to maintain the 
RMS pointing error within the required tolerance. 
The nominal structural mass was 442.04 kg and the 
actuator mass was assumed constant at 150 kg. 

With the average control power F ma x constrained 
at 3 N 2 -m 2 , the initial design gave an rms point- 
ing error of 73.6 /irad. The conventional design ap- 
proach was then followed, wherein the control gains 
(12 elements of the two Cholesky factors L p and 
L r ) were optimized for the fixed nominal structure. 
This control-optimized design yielded an rms point- 
ing error of 26.8 yurad. Next, an integrated design 


was performed, wherein both the structural and con- 
trol design variables were allowed to change simulta- 
neously. This integrated design resulted in an rms 
error of 16.8 /irad, which represented a 37-percent 
reduction over the conventional design. The struc- 
tural mass was slightly lower than the nominal de- 
sign. The upper and lower bound values, initial val- 
ues, and the optimal values of the structural design 
variables are summarized in table 3 and figure 5. The 
integrated design redistributed the structural mass 
from the battens and diagonals of the last two sec- 
tions of the main bus (closest to the small antenna) 
and small antenna support members to the large an- 
tenna support members and the section of the main 
bus closest to the large antenna, thus increasing the 
stiffness of these sections. This behavioral trend may 
be attributed to a trade-off between structural con- 
trollability and observability and its excitability by 
disturbances. In other words, the stiffness (or flexi- 
bility) of the structure is redistributed to establish a 
balance between the ability of the control system to 
fine-point the structure efficiently and the ability of 
the structure to reject disturbances. 

The modal frequencies of the first 10 flexible 
modes of the redesigned EPS structure (redesigned 
through the integrated design process) are presented 
in table 4. The data indicate that the first three 
frequencies (modes 1-3) associated with the small 
antenna support structure have been reduced sub- 
stantially, mainly because the pointing performance 
of the small antenna was not included in the per- 
formance metric. On the other hand, the frequency 
of the modes associated with the large antenna sup- 
port structure has increased considerably (for exam- 
ple, the first modal frequency increased from 0.58 Hz 
(table 1) to 2.75 Hz (table 4)), thus making the 
large antenna and support structure less sensitive to 
disturbances at the inputs or the outputs. At the 
same time, the first flexible frequency of the main 
truss (mode 7) decreased from about 6.6 Hz (ta- 
ble 1) to about 6.0 Hz (table 4), making the main 
truss more sensitive to disturbances but more con- 
trollable (observable). The elements of the 3x3 
lower triangular Cholesky factorization matrices L p 
and L r are given in table 5 for the nominal con- 
troller design, the control-optimized design, and the 
integrated controls-structures design. The attitude 
and rate gain matrices may be computed from the 
Cholesky matrices by using equation (14). 


To evaluate the effect of varying the actuator mass 
in the integrated design process, the actuator mass 
was allowed to vary by relating it to the infinity 
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norms of the gain matrices (a worst-case scenario); 
that is, 


mass and the inverse of P c is minimized. The design 
problem is to minimize 


A/act = a \u\ x = a (IGrL |y r U + IG,^ |y p y (23) 

in which a is a constant scalar taken to be 1.005 
and j |oo denotes the infinity norm of | |. It is noted 
that [ypl^ and \y r \^ were assumed fixed a priori to 
be 0.1 rad and 0.1 rad/sec, respectively. For this case, 
the actuator mass increased from 150 kg to 298.7 kg, 
whereas the rms pointing error and the structural 
mass were essentially unaffected, which is because 
the structure is rather stiff and is not affected by 
small masses. 

Design Problem II 

The system configuration for the second design 
problem is illustrated in the block diagram shown 
in figure 6. Here, the attitude and rate sensors are 
assumed to be polluted with noise modeled as zero- 
mean, white-noise processes. Furthermore, no dis- 
turbances at the inputs are considered. With this 
configuration, the system equations become 


A/str T T^act 
Kt, + M act 


+ ( 1 -/ 3 ) 


1/Pc 

1/p? 


(27) 


with respect to structural and/or control design vari- 
ables, subject to a constraint on e 


t !™ ( Tr { E e (0 e r (0 }) < emax (28) 


as well as side constraints on structural member sizes 


^min ^ ^ ^max (29) 

Note that both the mass term and the term for the 
controlled performance in equation (27) are normal- 
ized with respect to their corresponding initial values. 
In equation (27), M str and M act denote the structural 
(truss) mass and the actuator mass, and the super- 
script o denotes the initial value of the corresponding 
variable. The covariance of e is computed from the 
steady-state covariance of the state vector; that is, 


z — Az + Bu 


(24) 


y p = C p z + S p v p ' 
y r = C^z -j- S r v r > 
yper = O = CperZ > 


(25) 


where v p and v r denote the attitude and rate sen- 
sors noise vectors, and e represents the rms pointing 
error vector at the large antenna. In this design prob- 
lem, a dual-objective formulation is pursued, wherein 
both the total mass (structure plus actuator) and 
a measure of the closed-loop performance are opti- 
mized. The closed-loop performance measure is cho- 
sen as the sum of the time constants of the closed- 
loop eigenvalues (the reciprocals of the absolute value 
of the real part of the closed-loop eigenvalues), which 
represents a measure of the transient performance 
of the system (i.c., the smaller the time constants 
are, the faster the system dissipates transient distur- 
bances). The closed-loop or controlled performance 
P c is the inverse of the sum of the time constants of 
the system, such that 


Pc = 



(26) 


Now, to minimize the total mass and maximize 
the closed-loop performance, an objective function 
J consisting of a convex combination of the total 


Pe = & ^ l E ( eeT )] = ^ ( C P erPz C per) (30) 

where P z = ^lim E(zz^) denotes the closed-loop 

steady-state covariance of the state that is deter- 
mined from the solution of the following Liapunov 
equation (ref. 25): 

A d P z + P z A^ - -BG p W p G^B T - BG r W r G^B r (31) 

in which the attitude measurement noise intensity 
matrix, W p — FJ(v p v^), and the rate measure- 
ment noise intensity matrix, W r = E(v r vJ) } are as- 
sumed to be diagonal (i.e., W p = I 3x3 (arcsec) 2 and 
W r — 2xI 3x3 (arcsec/sec) 2 ). The pointing tolerance 
e max in equation (28) is taken as 11 ^/rad. 

The parameter ff allows for a trade-off between 
the two objectives. As /? is varied from one to zero, 
more emphasis is placed on the closed-loop perfor- 
mance and less on the total mass or the cost. Further- 
more, to satisfy certain fine-pointing requirements, 
an upper bound constraint is placed on the rms error 
at the large antenna due to white measurement noise 
in the attitude and rate sensors. Additional side 
constraints are also placed on the structural design 
variables for safety and practicality concerns. Fig- 
ures 7(a)-(d) show, respectively, the behavior (nor- 
malized relative to the nominal design) of the objec- 
tive function J, controlled performance (a measure 
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of transient performance), structural mass, and ac- 
tuator mass as a function of the trade-off parame- 
ter j3. As (3 approaches zero, the design becomes 
performance or control dominated, as evidenced by 
the sharp increase in the controlled performance 
(fig. 7(b)) and the substantial increases in the struc- 
tural and actuator masses (figs. 7(c) and (d)). On 
the other hand, as (3 approaches one, the controlled 
performance diminishes considerably, and the struc- 
tural mass and actuator mass decrease, giving a mass 
or cost- dominated design. 

One may use these trade-off results to obtain 
a design that gives a good controlled performance 
with acceptable cost levels by choosing a right value 
for {3. For example, at /? = 0.50, table 6 gives 
the corresponding integrated design with a control- 
optimized (conventional) design. An initial design to 
achieve good rigid-body performance (same as design 
problem I) was first obtained. This design required 
an actuator mass of about 171 kg (computed from 
eq. (23)), which along with the nominal structural 
mass of 442.04 kg and payload mass of 5168 kg, re- 
sulted in a total spacecraft mass of 5781.04 kg. The 
controlled performance measure P c was about 0.003. 
Then, a control-optimized design was performed by 
optimizing the control variables only. This design 
resulted in a 44-percent increase in the control per- 
formance, but also increased the total mass by 1 per- 
cent (57.17 kg). An integrated design was next car- 
ried out, which resulted in a 234-percent increase 
in the control performance along with lower total 
mass over the initial design (6 percent less than 
the total mass of the control-optimized design, or 
about 353 kg). Integrated design reduced the struc- 
tural mass by 53 percent and decreased the actuator 
mass by 24 percent. 

The optimization data and results for the struc- 
tural design variables are given in table 7 and fig- 
ure 8. It is observed that the structural mass was 
reduced and redistributed (i.e., mass was taken from 
the main bus structure and was added to the an- 
tenna support members). In other words, most of 
the longerons, battens, and diagonals of the main bus 
structure were reduced to lower bound values, while 
the antenna support members increased drastically 
in size. The open-loop modal frequencies of the re- 
designed EPS structure (obtained through integrated 
design process) are presented in table 8, where it is 
seen that the first flexible frequency corresponding 
to the large antenna support structure increased to 
near 1.8 Hz, with a reduction of truss flexible frequen- 
cies by as much as 50 percent at the higher end of 
the spectrum. In this configuration, the antennas are 
less susceptible to disturbances at the inputs or the 


outputs, whereas the overall structure has become 
more controllable (observable). 

The closed-loop eigenvalues are presented in 
figures 9 11 for the nominal design, the control- 
optimized design, and the integrated controls- 
structures design, respectively. These figures indi- 
cate that integrated design is more effective than the 
conventional design in moving the closed-loop eigen- 
values of the system farther into the left-half plane 
(such that the system can dissipate transient distur- 
bances faster), thus clearly demonstrating the advan- 
tage of the integrated design. Also, the Cholesky fac- 
torization matrices and L r (associated with the 
attitude and rate gain matrices) are given in table 9 
for the nominal design, the control-optimized design, 
and the integrated design. 

Summarizing the results obtained from design 
problems I and II, both clearly show the advantage of 
integrated design over the conventional approach in 
that the integrated design produces a better overall 
design. The main advantage of integrated design 
is its ability to obtain a better overall design and 
not necessarily just a reduction of the total mass. 
Moreover, a comparison of the behavioral trends of 
the structural design variables for the first and second 
design problems indicates that such trends could vary 
considerably depending on the design objectives and 
constraints of the integrated design problem. 

Integrated Design of the CSI 
Evolutionary Model 

This section considers the application of the inte- 
grated design methodology to the CSI Evolutionary 
Model. The discussion on the application is in En- 
glish units to be consistent with the formal model of 
the structure (ref. 14). The Phase- Zero Evolution- 
ary Model, shown in figure 12, basically consists of a 
62-bay central truss (each bay 10 in. long), two ver- 
tical towers, and two horizontal booms. Using two 
cables as shown, the structure is suspended from the 
ceiling (about 840 in. above the main truss). A laser 
source is mounted at the top of one of the towers, and 
a reflector with a mirrored surface is mounted on the 
other tower. The laser beam is reflected by the mir- 
rored surface onto a detector surface 660 in. above 
the reflector. Eight proportional bidirectional gas 
thrusters, with maximum output force of 4.4 lb each, 
provide the input actuation, and collocated servo ac- 
celerometers provide output measurements. The fi- 
nite element model of the system has 3216 degrees 
of freedom; therefore, the bulk of the computational 
effort is required for the solution of the structural 
eigenvalue problem of that size. 
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The design model consisted of the first 20 modes 
of the structure, including 6 suspension modes (i.e., 
modes due to the suspension of the structure and 
gravity) and 14 flexible modes. A modal damping 
ratio of 0.1 percent was assumed. The first 20 modes 
of the nominal Phase-0 Model are presented in ta- 
ble 10. The first six modes that range from 0.12 Hz 
to 1.16 Hz are the suspension modes (or non-zero- 
frcquency-rigid-body modes). The first two bending 
modes (lateral and vertical) and the first torsional 
modes of the structure are modes 7, 8, and 9. The 
system configuration for this design problem is il- 
lustrated in block diagram shown in figure 13. Here, 
white-noise disturbances of unit intensity are applied 
to the structure at the same locations as the con- 
trol inputs, and no disturbances at the outputs are 
considered. With no appreciable sensor and actua- 
tor dynamics considered, the system equations are 
written as 


z — Az -f Bu T Bw 

(32) 

y r = C r z 1 

(33) 

Yper = O = C p Z J 


where the performance vector y per is an 8 x 1 vec- 
tor corresponding to the displacements at the eight 
accelerometer locations. Here, only velocity measure- 
ments are used for feedback, and they are obtained 
by integrating the accelerometer signals with the aid 
of wash-out filters. It is noted that position feedback 
is not necessary for asymptotic stability because the 
structure is open-loop stable. 

To perform the integrated design, the structure 
was divided into seven sections (see fig. 12), three 
sections in the main bus, and one section each for 
the two horizontal booms and the two vertical towers. 
Three structural design variables were used in each 
section, namely, outer diameters of the longerons, 
the battens, and the diagonals, making a total of 
21 structural design variables. The design was posed 
in the form of a nonlinear optimization problem, 
wherein the steady-state average control power in the 
presence of a white-noise input disturbance with unit 
intensity (i.e., Standard deviation intensity = 1 lbf) 
is minimized, with the steady-state rms position error 
vector at the eight accelerometer locations and total 
mass constrained to required values. That is, the 
problem solved was: 

Minimize 

Hm Tv[£{ u (*) u r (f)}] (34) 


subject to 

{ E [ e (0 e7 w]}) 3 < Gnax (35) 
and a constraint on the total mass 

Mtot — Tfmax (36) 

Here, e max was taken to be 1.7 in. and M max was 
chosen equal to the total mass of the nominal Phase-0 
structure at 2.0 lb-sec 2 /in. 

Static Dissipative Controller 

As mentioned earlier, only velocity measurements 
are used for feedback, which are obtained by pro- 
cessing the accelerometer outputs. Thus, the static 
dissipative controller is given by 

u — G r y r (37) 

Here, the rate gain matrix was an 8 x 8 diagonal 
matrix (i.e., eight control design variables, along with 
the 21 structural design variables, resulted in a total 
of 29 design variables for the design optimization). 

The results of the design optimization are sum- 
marized in table 11. A conventional or control- 
optimized design was first performed (with the struc- 
tural design fixed at the nominal values), which 
required an average steady-state control power of 
19.34 lb 2 . 

Next, an integrated design was performed, where- 
in the average control power was minimized with re- 
spect to both control and structural design variables. 
The results (table 11) indicate an average control 
power of 11.41 lb 2 , a reduction of about 40 percent 
in the average control power over the conventional 
design. The initial and final values of the structural 
design variables, along with the corresponding lower 
bound and upper bound values, are presented in ta- 
ble 12 and figure 14. Keeping in mind that the tube 
diameters of the nominal CSI Evolutionary Model 
structure are 0.367 in. for the longerons and battens 
and 0.349 in. for the diagonal, it is observed from ta- 
ble 12 and figure 14 that sections 1 and 2 of the main 
bus (farthest away from the reflector) and the laser 
tower are considerably stiffened, whereas the hori- 
zontal booms and the reflector tower became more 
flexible, partly to satisfy the mass constraint. Gener- 
ally, in those sections that showed an increase in stiff- 
ness, the longerons increased in size more than the 
diagonals and the battens, because they were most 
effective in changing the stiffness of a section. 
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The trends in table 12 may be attributed to a 
trade-off between structural controllability, observ- 
ability, and excitability. The areas near the sources 
of disturbance (actuator locations) were stiffened in 
order to reduce the sensitivity of the structure to ex- 
ternal disturbances at those locations, while ensur- 
ing that no appreciable loss of controllability and/or 
observability occurred (as observed from the tran- 
sient controlled performance values in table 11). The 
small reduction in the transient controlled perfor- 
mance of the integrated design over the conventional 
design in table 11 is mainly because no measure of 
the transient performance was included in the op- 
timization process. Moreover, although the tran- 
sient controlled performance was lower in the inte- 
grated design case, it is still quite acceptable. The 
modal frequencies of the first 20 modes of the re- 
designed Phase-0 structure (redesigned through the 
integrated design process) are presented in table 13. 
These frequencies indicate that the first six frequen- 
cies associated with the suspended structure have 
not been changed significantly, mainly because the 
changes in the structure can affect these frequencies 
only through changing the location of the center of 
mass of the structure and not directly as is the case 
for the flexible modes. On the other hand, the fre- 
quencies of the flexible modes, particularly the sec- 
ond and third flexible modes, have increased consid- 
erably (as much as 40 percent). The second flexible 
mode frequency increased from 1.68 Hz to 2.33 Hz 
(mode 8 in table 13) and the third flexible mode 
frequency from 2.08 Hz to 2.65 Hz (mode 9 in ta- 
ble 13), making these modes and the structure less 
sensitive to disturbances at the inputs. The diagonal 
elements of the rate gain matrix are given in table 14 
for the control-optimized design and the integrated 
controls-structures design. 

Dynamic Dissipative Controller 


The next two designs that were performed for the 
CSI Evolutionary Model were conventional (control- 
optimized) and integrated designs with a dynamic 
dissipative controller. Because the system has no 
zero-frequency modes, Gp and G r were taken to be 
zero. The dynamic dissipative controller represented 
by equations (10) and (11) was used with block- 
diagonal compensator state matrix A c (consisting 
of eight second-order blocks) and compensator input 
influence matrix B c as follows: 


" A c i 0 

o a c2 


L o o 


0 1 
0 


A c8 J 


(38) 
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0 

B c i 

... 0 

(39) 

B c = 
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0 

b c8 . 


where A cl and B< 

•yi (* 

1,2,.. 

. , 8) are, respectively, 


2x2 matrices and 2x1 vectors, defined as 


— 





0 

1 



(40) 


Furthermore, the weighting matrix Q in equation (12) 
is assumed to be diagonal; that is, 


Q = Diag (<7i, <72 5 • - * > <Zie) ( 41 ) 

Here, the scalar variables (i = 1,2,..., 8), 

and qj (j = 1 , 2 ,..., 16) were chosen for the con- 
trol design variables. Thus, the number of control 
design variables was 32, making the total number of 
design variables 53. Table 15 shows the results of 
the designs with the dynamic dissipative controller. 
The conventional design reduces the control power 
by over 50 percent (to 8.31 lb 2 ), compared with the 
static dissipative case. The integrated design reduces 
the average control power by approximately another 
30 percent (over the conventional design) to 5.91 lb 2 . 
Table 15 also indicates that the transient controlled 
performance is slightly lower for the integrated design 
case. Similar to the static dissipative design, this re- 
duction in performance is because the optimization 
process did not include any measure of transient per- 
formance. Moreover, although tables 11 and 15 show 
that the transient controlled performance of the dy- 
namic dissipative design was quite lower than the 
static dissipative design, one should keep in mind 
that in the dynamic dissipative design, this perfor- 
mance measure includes contributions from the com- 
pensator poles, and thus it is not readily comparable 
with the static dissipative values. 

The optimization data and results for the struc- 
tural design variables are summarized in table 16 and 
figure 15. The trends are quite similar to the static 
dissipative design given in table 12. In fact, the two 
optimal structures are within 20 percent of each other 
(i.e., the structure that is optimal for the static dis- 
sipative design is also optimal for the dynamic dis- 
sipative design). The modal frequencies of the opti- 
mal structure for the dynamic dissipative design are 
summarized in table 17. These frequencies are also 
close (within 5 percent) to their corresponding values 
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for the static dissipative design. Finally the singu- 
lar value plots for the input-to-output transfer func- 
tions of the dynamic dissipative compensator are pre- 
sented in figures 16(a) -(h) for the control-optimized 
design and the integrated controls-structures design. 
The singular value plots for the control-optimized de- 
sign and the integrated design are somewhat similar 
except for input/output channel nos. 1 and 7 (see 
figs. 16(a) and (g)), where considerable reduction in 
the singular values, mainly in the lower end of the 
frequency spectrum, and increase in the bandwidth 
arc observed in the integrated design case over the 
conventional design. The elements of the compen- 
sator state matrix A c arc given in table 18 for the 
control-optimized design and the integrated controls- 
structures design. Also, the elements of the weighting 
matrix Q are presented in table 19 for both designs. 

The results obtained for both the static and 
dynamic dissipative controllers clearly show that 
integrated controls-structures design methodology 
can yield an overall design that is substantially supe- 
rior than the conventional sequential design scenario. 

Concluding Remarks 

An optimization-based approach has been devel- 
oped for performing integrated controls-structures 
design of a class of flexible spacecraft. The approach 
formulates the problem as a constrained optimization 
problem, wherein the set of design variables consists 
of both control and structural design variables. The 
approach uses static and dynamic dissipative control 
laws, which provide robust stability in the presence 
of parametric and nonparametric uncertainties in the 


model. The approach was demonstrated by appli- 
cation to integrated designs of a generic model of 
a space platform concept, as well as to a model of 
the Controls-Structurcs Interaction (CSI) program’s 
Evolutionary Model, which is a ground-based experi- 
mental testbed developed and constructed at Langley 
Research Center. 

The numerical results obtained indicate that the 
integrated design approach can yield substantially 
superior spacecraft design as compared with the tra- 
ditional sequential design approach. For example, a 
40-percent reduction in the pointing error is observed 
along with a slight reduction in mass, or an almost 
twofold increase in the controlled performance is in- 
dicated with more than a 5-percent reduction in the 
overall mass of the spacecraft (a reduction of hun- 
dreds of kilograms). Furthermore, the automated 
nature of the integrated design approach can accom- 
modate a wide variety of design specifications and re- 
quirements. A practical software tool (CSI-DESIGN) 
is being developed for performing integrated designs. 
Research is presently in progress for incorporating 
other types of control methods (such as or Tfocj 
into the integrated design process. Future plans also 
include hardware validation of the integrated design 
approach by constructing a laboratory test article 
based on the redesigned CSI Evolutionary Model 
and comparing the performance of the redesigned 
structure with the original CSI Evolutionary Model. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
October 14, 1992 


12 



Appendix 

Dissipative Controllers 

This appendix briefly presents the main results 
concerning dissipative controllers for flexible struc- 
tures. Dissipative controllers use the passivity prop- 
erties of the system that result from the use of col- 
located and compatible actuators and sensors, such 
as torque actuators and attitude and rate sensors or 
force actuators and position and velocity sensors. In 
general, dissipative controllers are designed to con- 
trol both rigid and elastic modes. In certain cases, 
the zero- frequency rigid modes may not be present. 
This may occur, for example, when paired torquers 
are used in a balanced configuration to control only 
the elastic motion (ref. 4) or in the case of ground- 
based test articles. In such cases, no position feed- 
back is required; only velocity feedback is sufficient to 
accomplish robust stability, and the dissipative com- 
pensator degenerates to a positivity controller. In 
the following, both the static and dynamic dissipa- 
tive controllers are discussed. The concepts and the 
underlying theory are based largely on the original 
work of Popov (ref. 26). 

Static Dissipative Controllers 

The constant-gain or static dissipative control law 
is given by: 

u = -G p y p - G r y r (Al) 

where G p and G r are m x m symmetric, positive- 
definite, position and rate gain matrices, and y p and 
y r represent the rax 1 position and rate measurement 
vectors. This control law is known to give guaranteed 
asymptotic stability of the entire system consisting of 
both rigid and flexible modes regardless of unmod- 
eled elastic modes or parameter uncertainties. (It 
was also shown in ref. 4 that stability is maintained 
even if small imprecision exists in the collocation 
of the actuators and sensors.) Furthermore, if Gp 
and G r are diagonal, this robust stability property 
is carried over in the presence of: 

(1) Monotonically increasing actuator nonlineari- 
ties, rate sensor nonlinearities belonging to the 
[0, oc) sector, and position sensor nonlineari- 
ties belonging to the (0,oo) sector (a single- 
valued function 'ip(v) is said to belong to the 
(0, oc) sector if ip(Q) = 0 and u'ip(iz) > 0 for 
v ^ 0; ip is said to belong to the [0, oo) sector 
if vip{v) > 0), and 

(2) Stable actuator dynamics g a (s) = k/(s + a), 
provided that g v jg r < a, where g p and g r de- 
note the appropriate diagonal elements of G p 


and G r . A possible drawback of these con- 
trollers is that the performance can be lim- 
ited because of the structure of the controller. 
The matrices Gp and G r can be designed 
to minimize a quadratic performance function 
or to obtain closed-loop eigenvalues close to 
the desired locations in the least-square sense 
(ref. 4). 

When the zero-frequency modes are absent, G p 
can be zero, and G r = G^ > 0 is sufficient for 
asymptotic stability. Furthermore, the closed-loop 
system is robust to [0, oc)-sector actuator and sen- 
sor nonlinearities and first-order actuator dynamics 
(ref. 4). 

Dynamic Dissipative Controllers 

To obtain better performance while still retain- 
ing the guaranteed robustness to unmodeled dynam- 
ics and parameter uncertainties, a class of dynamic 
dissipative controllers was considered. The case in 
which both zero- frequency rigid modes and flexible 
modes are present is considered first. 

Assuming for simplicity that the plant has three 
(one per axis) torque actuators and collocated atti- 
tude and rate sensors, the 3x3 transfer function from 
the torque input to the attitude rate output is given 
by: 



s 2 + 2 pibJiS + uf 


(A2) 


where J is the moment-of-inertia matrix and 
and u>i denote the rotational mode shape vector, 
damping ratio, and natural frequency of the ith 
structural mode. The transfer function G'(s) is 
positive real (PR), as defined as follows (ref. 27): 


Definition 


A rational matrix-valued function T(s ) of the 
complex variable s is said to be PR if T(s) is real 
when s is real, and T(s) +T t (s*) > 0 for Re(s) > 0, 
where * denotes the complex conjugate. 

Scalar PR functions have a relative degree (i.e., 
the difference between the degrees of the denomi- 
nator and numerator polynomials) of —1, 0, or 1 
(ref. 28). It can also be shown that PR matrices 
have no transmission zeros or poles in the open right- 
half of the complex plane and that the poles on the 
imaginary axis are simple and have nonnegative def- 
inite residues. By applying the maximum modu- 
lus theorem, it can be shown that it is sufficient to 
check the positive semidefiniteness of T(s) only on 
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the imaginary axis (s = ju;,0 < w < oc). Suppose 
(A,B,C, D) is an nth-order minimal realization of 
T(s). Prom reference 27, a necessary and sufficient 
condition for T(s) to be PR is that there exists an 
n x n symmetric positive-definite matrix P and ma- 
trices W and L such that 


A t V + PA = -LL t 
C = B t P + W r L 
W T W - D + D r 


The input to the plant is given by: 
u = -VK 

Define 


A z = 


A K 0 

c K o 


> 

(A3) 

B z = 

B A 

. Dr 


J 


C z = [0 

13X3] J 


(A7) 


(AS) 


This result is generally known in the literature 
as the Kalman-Yacubovich lemma. A stronger con- 
cept along these lines is strictly positive-real (SPR) 
systems. However, there are several nonequivalent 
definitions of SPR, all of which require the system to 
have all poles in the open left half plane (ref. 29). For 
the purpose of this paper, we define a less restrictive 
class of strongly PR systems as follows: 

Definition 


We present the following stability result that was 
proved in reference 16. 

Theorem 1 

Suppose 

(i) The matrix Ax is strictly Hurwitz 

(ii) There exists an (nx + 3) x (n K + 3) matrix 
P z = P^ > 0 such that 


A rational matrix-valued function T(s) of the 
complex variable s is said to be strongly PR if T(s) 
is real when s is real, and T(s) + T T (s*) > 0 for 
Re(s) > 0. 

The obvious difference between this definition and 
the definition of PR systems is that > has been 
replaced by strict inequality. The difference between 
the strongly PR and SPR systems is that the latter 
have poles only in the open left half plane, whereas 
the former can have poles on the jc^-axis. 

The transfer function from u to is given by: 


G(s) = 


G'(s) 


(A4) 


A^P Z +P Z A Z — — Q z — — diag ,03x3) (A9) 

where is a 3xn matrix such that (L/f, Ax) 
is observable, and Lx(sl — A#) -1 B/£ has no 
transmission zeros in Re(s) > 0 

(iii) The matrix C z is defined as 

C z = Bjp z (A10) 

(iv) The equation K(s) = Cx(sl - Ax)~ l Bx + 
Dx has no transmission zeros at the origin. 

Under these conditions, the controller K(s) 
stabilizes G(s). 


It can be seen that G(s) is not PR. However, G r (s), 
the transfer function from u to y r , is PR. 

Definition 

The compensator K(s) is said to stabilize a plant 
G(s ) if the closed-loop system consisting of stabiliz- 
able and detectable realizations of K and G in the 
standard feedback configuration is asymptotically 
stable (a.s.). 

Suppose a controller K(s) is represented by the 
n/^th-order minimal realization: 


*K = A K x K + B K y p 

(A5) 

y K = C A -x a - + D k y p 

(A6) 


Remark 1. In theorem 1, if is an nx x nx 
nonsingular matrix, then the observability and 
minimum-phase properties in condition (ii) are 
satisfied and the closed-loop system is a.s. 

Remark 2. The controller K stabilizes the full 
plant (i.e., the system consisting of the rigid modes, 
the elastic modes, and the compensator state vector 
(x#) is a.s.). The asymptotic stability is guaranteed 
regardless of the number of modes in the model or 
parameter uncertainties. The order of K can be 
chosen to be any number >3. In other words, this 
result enables the design of a controller of any desired 
order, which robustly stabilizes G. A procedure 
for designing K is to choose Q z = diag (Qas 03x3), 
where Q# = > 0, and to choose a stable Ax 

and matrices B# and Cx so that equations (A9) 
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and (A10) are satisfied. Using equation (A8) and 
defining 


Pzl Pz2 



where P z i is an n x x n k matrix and P Z 3 is a 3 x 3 
matrix, conditions (ii) and (iii) of theorem 1 can be 
expanded as: 


P zl ^a' + A/f P zl + P z2 + ^A' P z2 — ^ 


P W + p z3 c a: = 0 

B A' P zl + ^ = 0 


(All) 


b£p z2 + d£p z3 = i 


In addition, P z must be positive definite. Because of 
the large number of free parameters (i.e., A#, B#, 
Ca, L#) it is generally not straightforward to 
use equation (All) to obtain the compensator, and 
this problem remains an area of continuing research. 
Another method is to use the following s-domain 
equivalent of theorem 1. (See reference 16.) 


Theorem 2 


The compensator K(s) stabilizes G(s) if K(s ) 
has no transmission zeros at s = 0 and K(s)/s is 
strongly PR. 


populated K(s), however, there appear to be no 
straightforward methods. 

The controller K(s ) (eqs. (A5) and (A6)) is not 
strictly proper because of the direct transmission 
term D f c . From a practical viewpoint, it is sometimes 
desirable to have a strictly proper controller because 
it attenuates sensor noise as well as high-frequency 
disturbances. Furthermore, the most common types 
of controllers, which include the Linear-Quadratic- 
Gaussian (LQG) as well as the observer-pole place- 
ment controllers, are strictly proper (i.e., they have 
a first-order roll-off). In addition, the realization in 
equations (A5) and (A6) does not use the rate mea- 
surement y r . The following result from reference 16 
states that K can be realized as a strictly proper 
controller wherein both y p and y r are used. 

Theorem 3 

The plant G{s) is stabilized by the controller K f 
given by 


*k = A K -x-K + [Bat - A k L L] 


y P 

Yr 


(A15) 


= C K X K (A16) 

where L is a solution of 


D K -C K h = 0 (A17) 


The condition that K(s)/s be strongly PR 
is sometimes much easier to check than the 
conditions of theorem 1. For example, let 
K(s) = diag[A' 1 (s),A' 2 (s), AT 3 (s)], where 


Ki{s) 


S 2 + + (hi 

1 s 2 + Q] jS + a Qj 


(A12) 


It is straightforward to show that K(s)/s is strongly 
PR if (for i = 1,2,3) fcj, <J0n ay, fa, and 0u are 
positive, and 

«1i - Pli > o (A13) 

aufloi ~ aoiPli > 0 (A14) 

For higher order AVs, the conditions on the poly- 
nomial coefficients are harder to obtain. One sys- 
tematic procedure for obtaining such conditions for 
higher order controllers is the application of Sturm’s 
theorem (ref. 28). Symbolic manipulation codes can 
then be used to derive explicit inequalities. The con- 
troller design problem can be subsequently posed as 
a constrained optimization problem that minimizes 
a given performance function. For the case of fully 


There are many possible solutions for L. The 
solution that minimizes the Frobenius norm of L is: 

L = C T K (c K C T K y 1 D k (A18) 

For the case with no zero-frequency modes, the dy- 
namic dissipative controller degenerates to a positiv- 
ity controller that uses the feedback of only the veloc- 
ity. For this case, the robust stability is guaranteed 
with an SPR dynamic compensator (refs. 30 and 31). 
In terms of the state equations, it is sufficient for 
robust stability that the compensator satisfies the 
Kalman-Yacubovich condition (iii) with (A, L) ob- 
servable. A simpler sufficient condition, which also 
involves fewer parameters, is that D x = 0 and L is 
a (lower or upper) triangular matrix. 

The problem of designing dissipative controllers 
that also provide optimal performance has been a 
subject of active research. It has been shown that 
the static dissipative controller minimizes a certain 
quadratic performance index (ref. 4). Methods based 
on eigenstructure assignment with dissipativity con- 
straints have also been recently proposed (ref. 32). 
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For the case with no zero-frequency modes, it was 
shown in reference 33 that the selection of LQ reg- 
ulator and estimator weighting matrices in a certain 
manner results in an LQG controller that is also SPR, 
and therefore robustly stabilizes the plant. For the 
case where zero- frequency modes are present, de- 
sign methods based on numerical minimization of 
the “distance” between a nominal LQG controller 
and dynamic dissipative controllers were suggested 
in reference 16. 


Another approach for the case with zero- frequency 
modes is a two- level controller wherein an inner loop 
static dissipative controller with position and veloc- 
ity feedback is used, and a dynamic dissipative (SPR) 
velocity feedback controller is used in the outer loop. 
This approach was discussed in reference 5. 

In spite of some progress, the problem of designing 
optimal dissipative controllers remains an area of 
continuing research. 


16 



References 

1. Balas, Mark J.: Trends in Large Space Structure Con- 
trol Theory: Fondest Hopes, Wildest Dreams. IEEE 
Trans. Autom. Control , vol. AC-27, no. 3, June 1982, 
pp. 522-535. 

2. Ketner, G. L.: Survey of Historical Incidences With 

Controls- Structures Interaction and Recommended Tech- 
nology Improvements Needed To Put Hardware in Space. 
PNL-6846 (Contract DE-AC06-76R10 1830), Battelle 
Memorial Inst., Mar. 1989. 

3. Newsom, Jerry R.; Layman, W. E.; Waites, H. B.; and 

Hayduk, R. J.: The NASA Controls- Structures Inter- 

action Technology Program. NASA TM- 102752, 1990. 

4. Joshi, S. M.: Control of Large Flexible Space Structures. 
Volume 131 of Lecture Notes in Control and Information 
Sciences, M. Thoma and A. Wyner, eds., Springer- Verlag, 
1989. 

5. Joshi, S. M.; and Maghami, P. G.: Dissipative Com- 
pensators for Flexible Spacecraft Control. Proceedings of 
the 1990 American Control Conference , Volume 2 , IEEE 
Catalog No. 90CH2896-9, American Automatic Control 
Council, 1990, pp. 1955 1961. 

6. Singiresu, S. Rao; Pan, Tzong-Shii; and Venkayya, 
Vipperla B.: Robustness Improvement of Actively Con- 
trolled Structures Through Structural Modifications. 
AIAA J., vol. 28, Feb. 1990, pp. 353 -361. 

7. Balakrishnan, A. V.: Combined Structures- Controls Op- 
timization of Lattice Trusses. Comput. Methods Appl. 
Mech. & Eng., vol. 94, no. 1, 1992, pp. 131-152. 

8. Bossi, J. A.; Hunziker, K. S.; and Kraft, R. H.: Integrated 
Control/ Structure Design. NASA CR- 182020, 1990. 

9. Maghami, P. G.; Walz, J. E.; Joshi, S. M.; and Armstrong, 
E. S.: Integrated Controls-Structures Design Method- 
ology Development for a Class of Flexible Spacecraft. 
Third Air Force/NASA Symposium on Recent Advances in 
Multidisciplinary Analysis and Optimization — A Collec- 
tion of Technical Papers, Sept. 1990, pp. 1-6. 

10. Maghami, P. G.; Joshi, S. M.; and Lim, K. B.: Inte- 
grated Controls-Structures Design: A Practical Design 
Tool for Modern Spacecraft. Proceedings of the 1991 
American Control Conference, Volume 2, IEEE Catalog 
No. 91CH2939-7, American Automatic Control Council, 
1991, pp. 1465 1473. 

11. Maghami, P. G.; Joshi, S. M.; Elliot, K. B.; and 
Walz, J. E.: Integrated Design of the CSI Evolution- 
ary Structure: A Verification of the Design Methodol- 
ogy. Paper presented at the Fifth NASA/DOD Controls- 
Structures Interaction Conference (Lake Tahoe, Nevada), 
Mar. 3-5, 1992. 

12. Woodard, Stanley E.; Padula, Sharon L.; Graves, 
Philip C.; and James, Benjamin B.: An Optimization 
Method for Controlled Space Structures With Variable 
Actuator Mass. Proceedings of the Fourth NASA/DOD 


Controls/Structures Interaction Technology Conference, 
Andrew’ D. Swanson, compiler, WL-TR-91-3013, U.S. Air 
Force, Jan. 1991, pp. 411-428. 

13. Ramler, J.; and Durrett, R.: NASA’s Geostationary Com- 
munications Platform Program. A Collection of Tech- 
nical Papers — AIAA 10th Communication Satellite Sys- 
tems Conference , Mar. 1984, pp. 613 621. (Available as 
AIAA-84-0702.) 

14. Belvin, W. Keith; Elliot, Kenny E.; Bruner, Anne; 

Sulla, Jeff; and Bailey, Jim: The LaRC CSI Phase-0 
Evolutionary Model Testbed: Design and Experimen- 

tal Results. Proceedings of the Fourth NASA/DOD 
Controls /Structures Interaction Technology Conference, 
Andrew D. Swanson, compiler, WL-TR-91-3013, U.S. Air 
Force, Jan. 1991, pp. 594 613. 

15. Onoda, Junjiro; and Haftka, Raphael T.: Approach to 
Structure/Control Simultaneous Optimization of Large 
Flexible Spacecraft. AIAA J., vol. 25, no. 8, Aug. 1987, 
pp. 1133-1138. 

16. Joshi, S. M.; Maghami, P. G.; and Kelkar, A. G.: Dy- 
namic Dissipative Compensator Design for Large Space 
Structures. A Collection of Technical Papers, Volume 1 — 
AIAA Guidance, Navigation and Control Conference, 
Aug. 1991, pp. 467-477. (Available as AIAA-91-2650-CP.) 

17. Lefschetz, Solomon: Stability of Nonlinear Control Sys- 
tems . Academic Press Inc., 1965. 

18. Jamshidi, M.; and Herget, C. J., eds.: Computer-Aided 
Control Systems Engineering. North Holland, c.1985. 

19. Fenves, Steven J.; Perrone, Nicholas; Robinson, Arthur R.; 
and Schnobrich, William C., eds.: Numerical and Com- 
puter Methods in Structural Mechanics. Academic Press, 
Inc., 1973. 

20. Belvin, W. Keith; and Park, K. C.: Computer Im- 

plementation of Analysis and Optimization Procedures 
for Control-Structure Interaction Problems. A Col- 
lection of Technical Papers — AIAA Dynamics Special- 
ist Conference, Apr. 1990, pp. 32-41. (Available as 
AIAA-90-1194-CP.) 

21. Vanderplaats, G. N.: ADS — A Fortran Program for Auto- 
mated Design Synthesis- Version 1.10. NASA CR- 177985, 
1985. 

22. Vanderplaats, Garret N.: Numerical Optimization Tech- 
niques for Engineering Design - With Applications. 
McGraw-Hill, Inc., c.1984. 

23. George, Alan; and Liu, Joseph W.H.: Computer Solution 
of Large Sparse Positive Definite Systems. Prentice-Hall, 
Inc., c.1981. 

24. Kenny, Sean P.; Hou, Gene J.; and Belvin, W. Keith: 
Eigensensitivity in Integrated Design. Proceedings of 
the Fourth NASA/DOD Controls /Structures Interaction 
Technology Conference, Andrew D. Swanson, compiler, 
WL-TR-91-3013, U.S. Air Force, Jan. 1991, pp. 31-40. 

25. Stengel, Robert F.: Stochastic Optimal Control Theory 
and Application. John Wiley & Sons, Inc., c.1986. 


17 


26. Popov, V. M.: Hyper stability of Control Systems. Springer- 
Verlag, 1973. 

27. Anderson, B. D. O.: A System Theory Criterion for 

Positive Real Matrices. SIAM J. Control , vol. 5, no. 2, 
May 1967, pp. 171-182. 

28. Van Valkenberg, M. E.: Introduction to Modem Network 
Synthesis . John Wiley & Sons, Inc., c.1960. 

29. Lozano-Leal, Rogelio; and Joshi, Suresh M.: Strictly 

Positive Real Functions Revisited. IEEE Trans. Autom . 
Control, vol. 35, no. 11, Nov. 1990, pp. 1243 1245. 

30. Benhabib, R. J.; Iwens, R. P.; and Jackson, R. L.: Sta^ 
bility of Large Space Structure Control Systems Using 
Positivity Concepts. J. Guid. & Control, vol. 4, no. 5, 
Sept. Oct. 1981, pp. 487-494. 


31. McLaren, M. D.; and Slater, G. L.: Robust Multivariable 
Control of Large Space Structures Using Positivity. J. 
Guid., Control, & Dyn., vol. 10, no. 4, July-Aug. 1987, 
pp. 393 400. 

32. Maghami, P. G.; Gupta, S.; and Joshi, S. M.: De- 

sign of Dissipative Low-Authority Controllers Using an 
Eigensystem Assignment Technique. Paper presented at 
1992 American Control Conference (Chicago, Illinois), 
June 24-26, 1992. 

33. Lozano-Leal, R.; and Joshi, S. M.: On the Design of 
Dissipative LQG-Type Controllers. Recent Advances in 
Robust Control, Peter Dorato and Rama K. Yedavalli, 
eds., Inst, of Electrical and Electronics Engineers, Inc., 
c.1990, pp. 251-252. 



Table 1 . Open-Loop Modal Frequencies for 
the Nominal EPS Model 


Mode 

Frequency, Hz 

r 1 

0.58 

2 

0.73 

3 

0.91 

4 

2.40 

5 

2.99 

6 

3.18 

7 

6.65 

8 

7.36 

9 

16.31 

10 

17.73 


Table 2. Integrated Design of the EPS Model for Design Problem I 


[Static dissipative controller] 


Design 

RMS 

pointing, 

/irad 

Structural 

mass, 

k g 

Actuator 

mass, 

kg 

Control 

power, 

N 2 *m 2 

Initial design 

73.6 

442.06 

150 

2.98 

Control-optimized design 

26.78 

442.06 

150 

3.00 

Integrated design (without actuator mass) 

16.78 

437.34 

150 

3.00 

Integrated design (with actuator mass) 

17.01 

400.32 

298.73 

3.00 


Table 3. Optimization Data and Results for the Structural Design Variables 
of the EPS Model for Design Problem I 




Upper bound 

Lower bound 

Initial 

Final 

Design variables 

Section 

values, m 

values, m 

values, m 

values, m 

1 (longeron) 

1 

0.15 

0.01 

0.051 

0.107 

2 (batten) 

1 

.15 

.01 

.051 

.030 

3 (diagonal) 

i 

.15 

.01 

.051 

.025 

4 (longeron) 

2 

0.15 

0.01 

0.051 

0.066 1 

5 (batten) 

2 

.15 

.01 

.051 

.010 

6 (diagonal) 

2 

.15 

.01 

.051 

.010 

7 (longeron) 

3 

0.15 

0.01 

0.051 

0.066 

8 (batten) 

3 

.15 

.01 

.051 

.041 

9 (diagonal) 

3 

.15 

.01 

.051 

.058 

10 (support) 

Large antenna 

0.15 

0.01 

0.051 

0.149 

11 (support) 

Small antenna 

0.15 

0.01 

0.051 

0.010 
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Table 4. Modal Frequencies for the Redesigned 
EPS Model for Design Problem I 


Mode 

Frequency, Hz 

1 

0.17 

2 

.22 

3 

.23 

4 

2.75 

5 

3.11 

6 

3.30 

7 

5.98 

8 

6.62 

9 

8.58 

10 

14.13 


i 

i 


Table 5. Elements of Cholesky Matrices for the Attitude and Rate Gain Matrices 
of the EPS Model for Design Problem I 


Design 

variable 

Nominal 

design 

Control-optimized 

design 

Integrated 

design 

Lp(l> 1) 

120.0 

188.7 

254.3 

Lp (2, 1) 

0 

7.1 

1.2 

L p (3,l) 

0 

3.9 

-1.7 

L P (2>2) 

180.0 

179.2 

245.0 

L p (3,2) 

0 

7.4 

0.6 

L P (3,3) 

240.0 

189.9 

255.7 

^(l, 1) 

90.0 

212.5 

223.2 

M2, 1) 

0 

-17.4 

-2.9 

L r (3,l) 

0 

15.3 

8.3 

L r (2,2) 

150.0 

272.1 

327.9 

M3, 2) 

0 

-33.0 

-1.0 

L r (3,3) 

200.0 

293.8 

403.3 


Table 6. Integrated Design of the EPS Model for Design Problem II (Static Dissipative Controller) 

[rms <11 //rad] 


Design 

Normalized 

controlled 

performance 

Normalized 

structural 

mass 

Normalized 

actuator 

mass 

Normalized 

total 

mass 

Initial design 

1.0 

1.0 

1.0 

1.0 

Control-optimized design, (3 — 0.50 

1.44 

1.0 

1.33 

1.01 

Integrated design, /? = 0.50 

2.34 

.47 

.76 

.95 
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Table 7. Optimization Data and Results for the Structural Design Variables 
of the EPS Model for Design Problem II 




Upper bound 

Lower bound 

Initial 

Final 

Design variable 

Section 

values, m 

values, m 

values, m 

values, m 

1 (longeron) 

1 

0.15 

0.01 

0.051 

0.010 

2 (batten) 

1 

.15 

.01 

.051 

.010 

3 (diagonal) 

1 

.15 

.01 

.051 

.010 

4 (longeron) 

2 

0.15 

0.01 

0.051 

0.018 

5 (batten) 

2 

.15 

.01 

.051 

.010 

6 (diagonal) 

2 

.15 

.01 

.051 

.010 

7 (longeron) 

3 

0.15 

0.01 

0.051 

0.010 

8 (batten) 

3 

.15 

.01 

.051 

.010 

9 (diagonal) 

3 

.15 

.01 

.051 

.010 

10 (support) 

Large antenna 

0.15 

0.01 

0.051 

0.114 

11 (support) 

Small antenna 

0.15 

0.01 

0.051 

0.083 


Table 8. Modal Frequencies for the Redesigned 
EPS Model for Design Problem II 


Mode 

Frequency, Hz 

1 

1.79 

2 

2.04 

3 

2.16 

4 

3.48 

5 

3.76 

6 

4.63 

7 

5.52 

8 

5.74 

9 

9.26 

10 

9.71 


Table 9. Elements of Cholesky Matrices for the Attitude and Rate Gain 
Matrices of the EPS Model for Design Problem II 


Design 

variable 

Nominal 

design 

Control-optimized 

design 

Integrated 

design 

MM) 

120.0 

201.9 

178.4 

M 2 , i) 

0 

1.2 

-.005 

Lp(3, 1) 

0 

1.5 

-.001 

M2, 2) 

180.0 

254.3 

178.4 

M 3 , 2) 

0 

1.8 

.012 

ljp(3, 3) 

240.0 

254.6 

178.3 

MM) 

90.0 

184.0 

206.0 

M2,i) 

0 

4.3 

.001 

M3, 1) 

0 

-1.8 

.002 

M2, 2) 

150.0 

251.5 

205.5 

M3, 2) 

0 

-2.4 

-.006 

M3, 3) 

200.0 

252.0 

200.9 
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Table 10. Modal Frequencies for the 
CSI Phase-0 Evolutionary Model 


Mode 

Frequency, Hz 

i 

1 

0.122 


2 

.126 

I 

3 

.173 

i 

4 

.680 

i 

5 

.704 

* 

6 

1.159 

{ 

7 

1.572 


8 

1.676 


9 

2.085 


10 

3.867 

= 

11 

3.972 

i 

12 

4.127 

_ 

13 

4.172 

i 

14 

5.788 

1 

15 

6.456 

| 

16 

6.568 


17 

6.777 

| 

18 

7.806 

1 

19 

8.732 


20 

9.396 



Table 11* Integrated Design of the CSI Evolutionary Model 
[Static dissipative controller] 



rms 

Control 

Transient controlled 

Design 

displacement, in. 

power, lb 2 

performance, P c 

Control-optimized design 

1.70 

19.34 

0.0048 

Integrated design 

1.70 

11.41 

.0038 
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Table 12. Optimization Data and Results for the Structural Design Variables of the CSI Evolutionary Model 

[Static dissipative controller] 




Upper bound 

Lower bound 

Initial 

Final 

Design variable 

Section 

values , in. 

values, in. 

values, in. 

values, in. 

1 (longeron) 

1 

1.0 

0.12 

0.3125 

0.999 

2 (batten) 

1 

1.0 

.12 

.3125 

.633 

3 (diagonal) 

i 

1.0 

.12 

.3125 

.958 

4 (longeron) 

2 

1.0 

0.12 

0.3125 

0.124 

5 (batten) 

2 

1.0 

.12 

.3125 

.644 

6 (diagonal) 

2 

1.0 

.12 

.3125 

.144 

7 (longeron) 

3 

1.0 

0.12 

0.3125 

0.999 

8 (batten) 

3 

1.0 

.12 

.3125 

.453 

9 (diagonal) 

3 

1.0 

.12 

.3125 

.658 

10 (longeron) 

4 

1.0 

0.12 

0.3125 

1.000 

11 (batten) 

4 

1.0 

.12 

.3125 

.181 

12 (diagonal) 

4 

1.0 

.12 

.3125 

.407 

13 (longeron) 

5 

1.0 

0.12 

0.3125 

0.425 

14 (batten) 

5 

1.0 

.12 

.3125 

.143 

15 (diagonal) 

5 

1.0 

.12 

.3125 

.138 

16 (longeron) 

6 

1.0 

0.12 

0.3125 

0.122 

17 (batten) 

6 

1.0 

.12 

.3125 

.392 

18 (diagonal) 

6 

1.0 

.12 

.3125 

.145 

19 (longeron) 

7 

1.0 

0.12 

0.3125 

0.142 

20 (batten) 

7 

1.0 

.12 

.3125 

.183 

21 (diagonal) 

7 

1.0 

.12 

.3125 

.427 
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Table 13. Modal Frequencies for the Redesigned 
CSI Evolutionary Model 

[Static dissipative controller] 


Mode 

Frequency, Hz 

1 

0.122 

2 

.127 

3 

.173 

4 

.625 

5 

.684 

6 

1.191 

7 

1.642 

8 

2.334 

9 

2.651 

10 

3.829 

11 

3.995 

12 

4.226 

13 

4.598 

14 

5.379 

15 

7.092 

16 

7.639 

17 

7.910 

18 

8.599 

19 

8.994 

20 

9.020 


Table 14. Diagonal Elements of Rate Gain Matrix for the CSI Evolutionary Model 

[Static dissipative controller] 


Design 

variable 

Control-opt imized 
design 

Integrated 

design 

G r (l, 1) 

1.525 

1.263 

G r (2,2) 

.409 

.477 

G r (3,3) 

1.234 

1.203 

G r (4, 4) 

.448 

.440 

G r (5, 5) 

2.126 

1.674 

G r (6, 6) 

2.122 

1.671 

G r (7, 7) 

.352 

.869 

G r (8, 8) 

.680 

.626 





Table 15. Integrated Design of the CSI Evolutionary Model 
[Dynamic dissipative controller] 


Design 

RMS 

displacement, in. 

Control 
power, lb 2 

Transient controlled 
performance, P c 

Control-optimized design 
Integrated design 

1.70 

1.70 

8.31 

5.91 

0.0026 

.0024 


Table 16. Optimization Data and Results for Structural Design Variables of the CSI Evolutionary Model 

[Dynamic dissipative controller] 


Design variable 

Section 

Upper bound 
values, in. 

Lower bound 
values, in. 

Initial 
values, in. 

Final 
values, in. 

1 (longeron) 

2 (batten) 

3 (diagonal) 

1 

1 

1 

1.0 

1.0 

1.0 

0.12 

.12 

.12 

0.3125 

.3125 

.3125 

0.999 

.547 

1.000 

4 (longeron) 

5 (batten) 

6 (diagonal) 

2 

2 

2 

1.0 

1.0 

1.0 

0.12 ! 
.12 
.12 

0.3125 

.3125 

.3125 

0.120 

.540 

.120 

7 (longeron) 

8 (batten) 

9 (diagonal) 

3 

3 

3 

1.0 

1.0 

1.0 

0.12 

.12 

.12 

0.3125 

.3125 

.3125 

1.000 

.354 

.579 

10 (longeron) 

11 (batten) 

12 (diagonal) 

4 

4 

4 

1.0 

1.0 

1.0 

0.12 

.12 

.12 

0.3125 

.3125 

.3125 

1.000 

.138 

.367 

13 (longeron) 

14 (batten) 

15 (diagonal) 

5 

5 

5 

1.0 

1.0 

1.0 

j 0.12 

.12 
.12 

0.3125 

.3125 

.3125 

0.357 

.120 

.120 

16 (longeron) 

17 (batten) 

18 (diagonal) 

6 

6 

6 

1.0 

1.0 

1.0 

0.12 

.12 

.12 

0.3125 

.3125 

.3125 

0.120 

.308 

.120 

19 (longeron) 

20 (batten) 

21 (diagonal) 

7 

7 

7 

1.0 

1.0 

1.0 

0.12 

.12 

.12 

0.3125 

.3125 

.3125 

0.120 

.135 

.332 


25 




Table 17. Modal Frequencies for the Redesigned 
CSI Evolutionary Model 

[Dynamic dissipative controller] 


Mode 

Frequency, Hz 


1 

0.122 


2 

.127 


3 

.173 


4 

.629 


5 

.690 


6 

1.186 


7 

1.578 


8 

2.238 


9 

2.576 


10 

3.761 

i 

11 

3.782 

I 

12 

4.174 

! 

13 

4.320 

| 

14 

5.196 

i 

15 

6.443 

‘ 

16 

7.351 


17 

7.876 

= 

18 

8.549 

| 

19 

8.922 

j 

20 

8.982 

| 


Table 18. Control Design Variables for the CSI Evolutionary Model 
[Dynamic dissipative controller] 


Design 

variable 

Control-optimized 

design 

Integrated 

design 

A c (2, 1) 

11.400 

45.204 

A c (2, 2) 

66.351 

107.918 

A c (4, 3) 

82.667 

72.342 

Ac (4, 4) 

92.846 

80.691 

Ac (6, 5) 

20.989 

35.454 

A c (6, 6) 

70.983 

94.716 

A c (8, 7) 

89.514 

74.961 

A c (8, 8) 

92.700 

79.441 

A c (10, 9) 

27.749 

43.161 

A c (10, 10) 

42.645 

88.555 

A c (12, 11) 

21.109 

44.916 

A c (12, 12) 

50.887 

91.610 

A c (14, 13) 

21.495 

38.261 

A c (14, 14) 

100.329 

86.221 

A c (16, 15) 

22.814 

29.081 

A c (16, 16) 

81.224 

90.669 
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Table 19. Control Design Variables for the CSI Evolutionary Model 
[Dynamic dissipative controller] 






i 


Figure 2. Flow diagram for the integrated controls-structures optimization. 
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Figure 3. Schematic of the EPS. 



Figure 4. Block diagram of the system configuration for design problem I. 
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Figure 6. Block diagram of the system configuration for design problem II. 
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(a) Objective function versus /?. (b) Normalized controlled performance versus /?. 



P P 


Performance dominated 


Mass dominated 


■ Performance dominated 


Mass dominated 


(c) Normalized structural mass versus /?. 


(d) Normalized actuator mass versus 3 . 


Figure 7. Controls-structures trade-off study for the EPS Model. 
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Figure 8. Optimized structural parameters of the EPS Model for design problem II. 



Figure 9. Closed-loop eigenvalues of the EPS model for design problem II nominal controller design. 
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Figure 10. Closed-loop eigenvalues of the EPS model for design problem II — control-optimized design. 
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Figure 11. Closed-loop eigenvalues of the EPS model for design problem II — integrated controls-structures 
design. 
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Figure 14. Optimized structural parameters for the CSI Evolutionary Model — static dissipative controller. 


/ longeron 



Figure 15. Optimized structural parameters for the CSI Evolutionary Model — dynamic dissipative controller. 
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(b) Input /output channel 2. 


Figure 16. Singular value plots for input/output channels of the dynamic dissipative controller -the CSI 
Evolutionary Model. 
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(c) Input/output channel 3. 




-60 


Conventional 

Integrated 


-70 h 


80 


J ■ i i i ml 

.001 


I i i ' ' ' m i ' i'i 

.01 .1 

Frequency, Hz 


(d) Input/output channel 4. 
Figure 16. Continued. 
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(g) Input /output channel 7. 
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(h) Input/output channel 8. 
Figure 16. Concluded. 
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